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ABSTRACT 

We present an astrometry and photometry catalogue of globular cluster (GC) candidates detected 
with HST WFPC2 in a sample of 19 early-type galaxies, appropriate for comparison to the low- 
mass X-ray binary (LMXB) populations observed with Chandra. In a companion paper, we present 
the Chandra data and investigate the relation between these populations. We demonstrate that, 
although there is little evidence of a colour-magnitude correlation for the GCs, after estimating 
mass and metallicity from the photometry, under the assumption of a single age simple stellar pop- 
ulation, there is a significant positive correlation between mass and metallicity. We constrained 
[Z/H] = (-2.1 ± 0.2) + (0.25 ± 0.04)logi M, with a 1-a intrinsic scatter of 0.62 dex in metallicity. If 
GCs are bimodal in metallicity this relation is consistent with recent suggestions of a mass-metallicity 
relation only for metal-poor clusters. Adopting a new technique to fit the GC luminosity function 
(GCLF) accounting for incompleteness and the Eddington bias, we compute the V-band local GC 
specific frequency (Sn) and specific luminosity (Sl) of each galaxy. We show that Sl is the more 
robust measure of the richness of a GC population where a significant fraction is undetected due to 
source detection incompleteness. We find that the absolute magnitude of the GCLF turnover exhibits 
intrinsic scatter from galaxy to galaxy of ~0.3 mag (1-c), limiting its accuracy as a standard distance 
measure. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: star clusters — galaxies: general — 
galaxies: distances and redshifts 
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1. INTRODUCTION 

Globular clusters (GCs) are found in galaxies of all 
morphological types and sizes. As primarily old stellar 
systems, their distribution and properties provide crucial 
insights into the way in which structure forms within 
th e Universe. For a recen t review we refer the reader 
to Brodie & Strader (2006). An intriguing characteris- 
tic oF+Eese^Ejee^Sy^vEel*^ 



field form i n the field, despite the stellar population being 
old. Juett (2005) reported a similar result although, since 
they did not strictly compare the numbers of LMXBs and 
GCs within the same apertures, the str ength of their cor - 
relation has been called into question ( [Kim et aL 200C| ). 
Typically ~4% of GCs are observed to harbour LMXBs, 
with brighter and redder (mo re metal rich) GCs preferen- 
tially likely to contain them ( Kundu ct al.| 2002| ; Barazin 



thcirl structure and internal dynamics, is their associa gt_aL| |200j; |Kim ct al.||2006|). The GC luminosity de- 

pendence is broadly consistent with the probability of 



tion with low-mass X-ray binaries (LMXBs). Although 
it has long been recognized that LMXBs are overabun- 
dant ipeMM«-k-eptie&14ight-by-*fe^ 



harbouring an LMXB being proportional to the mass or 
to the stellar capture cross-section in the GC (e.g. Jordan 



GCslasVnmpRreHTn thefipld"(}Fabimi gj^Sgkl e jQS B EOsfe et g B jivakoff et gp| The 



1975), the small number of sources has limited what this 
phenomenon can tell us about LMXB formation. With 



origin of the colour dependence is unclear and several 
different possible explanations have been propos e d (for 



dvont of Chandra, however, it hao become 



the 

to rebolve indmduaJ LMXBs in galaxies nntsiHe the Ln- 



fc acib lc a review of s o me of these, see Jordan et al. 2004 ) 



Mac- 



cal Group, and thus to begi n to assemble larger samples 



of LMXB-hosting GCs (e.g. Barazin et al 



2003[ ). Given 

their typically rich GC populations and their clean, old 
stellar populations which prevent contamination of the 
X-ray sources with high-mass X-ray binaries, massive 
early-type galaxies, in particular, provide an i deal en 
vironment in which to conduct such studies (e.g . [Kundu 



carone et al. (2004) suggested it may relate to different 
patterns of mass-transfer in metal-rich and metal-poor 
LMXBs. Alternatively, it may arise due to variations 
i n the IMF betw een metal-poor and metal-rich clusters 



FFal 



J2002 
Even wit 



arazin ct al. 2003] ; Humphrey & Buotc 2004 f. 
Dies o 



the small samples of galaxies in which the 



GC-LMXB connection has been investigated to date, 
some intriguing trends have already been observed. By 
comparing the numbers of GCs an d the t otal l uminos- 
ity of LMXBs in early- type galaxies Irwin ( [2005 ) argued 
that a significant fraction of the LMXBs observed in the 

1 rWiart.mpnt nf Ptiygirs and Agtrnnnmy University nf Halifnr- 



nia, i jvine, 4129 Frederick Reines Hall, Irvine, CA 92697-4575 



( Grindlay 1987), or the efficiency of magnetic braking 
(Ivanova 2006 J! So far, however, most studies of the 
LMXB-GC connection have been relatively small and so 
strong general conclusions are difficult to draw. 

The properties of the GC populations in early-type 
galaxies themselves are also of considerable interest since 
they provide a unique insight into how such galaxies 
form. Although the numbers of GCs do not simply scale 
with the total stellar mass of a galaxy, the GC luminosity 
func tion (GCLF ) is observed to be remarkably uniform 
(e.g. Harris 1991 ). Typically it can be well-fitted by a log- 
normal distribution, the absolute shape of w hich varies 
only weakly with the galaxy proper ties (e.g 
Whitmorc||2001a| ; Jordan et alj f2006|) 



Kundu 
The origin of this 
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shape may be relat ed to the dynamical d estruction of 
low-mass GCs (e.g. Vesperini fc Zcpl 2003 ). The stabil- 
ity of the GCLF peak ( "turnover" ) has led to its ado ption 



rig 



as a standard distance measu re for nearby galaxies ( Har ■ 
1391; Jacoby et al. 1992j ), although there has been 



some debate as to its reliability (e.g. Ferrarese et al 
|2000| ; |Kundu fe Whitmorc| |2001a| ). A recent study of 
early-typ e galaxies in Virgo (which span ~7 magnitudes 
in B) by Jordan et al. (2007) found that the width and, 
possibly, the turnover magnitude vary with the galaxy 
absolute magnitude, in the sense that the least-massive 
galaxies have narrower, fainter GCLFs. The authors did, 
however, find significant scatter about these relations for 
galaxies at a given magnitude. 

The GC colour distributions of early-type gala xies are 
obse rved to be almost universally bimodal, (e.g . Zepf & 
Ashr fianl |l993| ; |Gcbhardt fe Kissler-Patigl |1999| ), proba- 



bly indicating ubiquitous me tal-rich and m e tal-poor sub 



Richtlci 


200G 




Yoon et al. 



[2006 ). The GC distribution, in particular that of the 



metal-poor (i.e. blue) population, is observed to be sub 



stantia lly m ore extended than the optical light (|Gcislcr 
et aQ |1996| ). The average colour of the GCs in a 
galax y strongly correlates wit h the total galaxy magni- 
tude ( [Brodie fe Huchra 1991), reflecting an increasing 
fraction of red GCs in more massive galaxies coupled with 
a correlation between mean GC met allicity and galax y 
mass for both blue and red clusters (Peng et al. 2006). 
There is, however, no correspondingly strong correlation 
betw een the colour and l uminosity of the GCs themse lves 



(e.g. Larsen et al. 2001; Kundu fc Whitmorc 2001a). A 



weak colour-magnitude relation has been observed, but 
only in m etal-poor GCs, in recent observations o f nearby 
galax ies ( | Harris et al. 2006; Strader et al. 2006 ; Mioskc 
et al. 2006 ) . These clues point to different origins for the 
two sub-populations, with the metal-poor clusters possi- 
bly forming in low-mass halos in the very early Universe 
and being accreted during hierarchical structure forma- 
tion. In contrast, metal rich clusters may form during 
t he gas-rich mergers wh ich assemble the massive galaxy 
( Brodie fc Strader 2006 , and references therein). 

In this paper, we present a study of the GC pop- 
ulations in 19 early- type galaxies observed with HST 
WFPC2. In particular, we present a photometric cat- 
alogue of GC candidates which can be directly compared 
to the properties of observ e d LM XBs. In a companion 
paper, Humphrey & Buotc (200S, hereafter Paper I), we 
present the Chandra data for these galaxies and investi- 
gate the relation between these populations, significantly 
expanding the sample of galaxies in which the LMXB-GC 
connection has been investigated. In the present work we 
compute self-consistently derived properties of the GC 
populations, such as the GC specific frequency, Sn, spe- 
cific luminosity, Sl, and the GCLF turnover, mj. The 
self-consistent computation of such parameters is essen- 
tial for our analysis in Paper I. The galaxies used in this 
study were chosen from the Chandra and HST archives 
as having sufficiently deep ACIS and WFPC2 data to 
enable a significant fraction of the LMXBs and GCs to 
be individually resolved. To enable a consistent com- 
parison between the galaxies and to mitigate against the 
rapidly-degrading Chandra PSF off-axis, we focus here 
only on WFPC2 pointings of the centre of each galaxy. 
The galaxies, and the details of the archival observations 



used are shown in Table [p. All errors quoted are 90% 
confidence regions, unless otherwise stated in the text. 

2. OBSERVATIONS AND DATA ANALYSIS 

Data reduction and analysis were performed using the 
Image Reduction and Analysis Facility, IRAF (PC-IRAF 
vers. 2.12.2a), and the Space Telescope Science Data 
Analysis System, STSDAS (vers. 3.2) software suites. 
The individual pipelined images were first processed to 
correct for warm pixels with the warmpix STSDAS task. 
The corrected images were then aligned and, where the 
gain and filters were identical, combined with the STS- 
DAS task gcombine to remove cosmic ray events. The 
alignment was achieved initially by examining the as- 
trometry keywords in the data-file and shifting individ- 
ual images an appropriate integer number of pixels. In 
a few cases the astrometry keywords were insufficiently 
accurate and so the offset between the images was de- 
termined by matching moderately bright point-sources 
between the exposures. In order to remove residual cos- 
mic ray events and defective pixels which have not been 
removed in earlier processing (in particular since there 
are some pointings for which only a single exposure is 
available), the resulting i mages were processe d with the 
LA-Cosmic 2 algorithm of van Dokkum (2001). 

The detection of point sources against the bright, 
strongly-varying, "background" of galaxy emission is 
particularly challenging. To achieve this, we adopted a 
procedure similar to that ou tlined by Kundu et al. (199E, 
Kundu & Whitmorc 2001a| ). First, to improve S/N, for 
each WFPC2 CCD we aligned and combined the images 
accumulated with different filters. An initial source list 
was created with the IRAF task daofind, having first es- 
timated the statistical noise level from the region of the 
CCD with the least background galaxy emission and free 
from obvious point-sources. Due to the rapidly rising 
background across the WFPC2 field of view this pro- 
cedure yields a very large number of detections, most 
of which are false on each CCD. To eliminate false de- 
tections, we recomputed the S/N for each source, based 
on local background estimates. To measure the source 
counts, we adopted a 3 pixel radius aperture and we es- 
timated the background from an annulus with inner and 
outer radii of 3 and 5 pixels. We rejected any source 
with S/N lower than 4.8, implying only ~1 false detec- 
tion per CCD. In order to provide useful colour informa- 
tion, we required that any source was detected in all fil- 
ters used. To eliminate spatially extended sources (most 
likely background galaxies), we excluded any source with 
a "concentration" (defined as the ratio of background- 
subtracted counts within a 3 pixel aperture to that within 
a 0.5 pixel aperture) outside the range 2-10. All detec- 
tions were confirmed by visual inspection of the images. 

To minimize confusion with small-scale but extended 
features (such as central dust-lanes) in the galax- 
ies, which seriously complicate the computation of 
background-subtracted photometry, we eliminated detec- 
tions within the central few arcseconds of each galaxy 
(shown in Table Q). In particular, data in the vicin- 
ity of a central nuclear disc in NGC 1387 and the bright 
knots of emission in the centre of NGC 4125 were ex- 
cluded. The projection of the bright spiral NGC 4647 

2 Available at http://www.astro.yale.edu/dokkum/lacosmic 
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TABLE 1 
Target list and observation log 



Name 


Type 


Dist 


L K 


B 




V 






I 
















Filt ObsID 


Exp 


Filt 


ObsID 


Exp 


Filt 


ObsID 


Exp 


(") 






(Mpc) 


(lO n L ) 




(s) 






(s) 






(s) 




NGC1332 


S(s)0 


21.3 


1.4 






F555W 


u2tv03 


160 


F814W 


u2tv03 


320 


2.9 


NGC1387 


SAB(s)0 


18.9 


0.78 






F606W 


u29r5h 


160 








8.7 


NGC1399 


cD 


18.5 


2.1 


F450W u34m02 


2600 








F814W 


u34m02 


1200 


0.31 


NGC1404 


El 


19.5 


1.5 


F450W u34m04 


2400 








F814W 


u34m04 


1260 


1.7 


NGC1553 


SA(rl)0 


17.2 


1.9 






F555W 


u2tv28 


160 


F814W 


u2tv28 


320 


2.0 


NGC3115 


SO 


9.00 


0.74 






F555W 


u2j20b 


860 


F814W 


u2j20b 


1236 


2.7 


NGC3585 


E7/S0 


18.6 


1.5 






F555W 


u3m719 


800 


F814W 


u3m719 


1800 


3.7 


NGC3607 


SA(s)0 


21.2 


1.5 






F555W 


u2tvl2 


160 


F814W 


u2tvl2 


320 


8.4 


NGC4125 


E6pec 


22.2 


1.8 






F555W 


u3m723 


1400 


F814W 


u3m723 


2100 


7.1 


NGC4261 


E2-3 


29.3 


2.2 






F547M 


u2i502 


800 


F791W 


u2i502 


800 


1.6 


NGC4365 


E3 


19.0 


1.6 






F555W 


u2bm07 


1000 


F814W 


u2bm07 


460 


3.5 


NGC4472 


E2/S0(2) 


15.1 


3.2 






F555W 


u66306 


7200 


F814W 


u66315 


5200 


2.9 


NGC4494 


El-2 


15.8 


0.81 






F555W 


u3vz05 


1500 


F814W 


u3vz05 


1800 


2.3 


NGC4552 


E 


14.3 


0.85 






F555W 


u30708 


1800 


F814W 


u30708 


1500 


0.85 


NGC4621 


E5 


17.0 


1.2 






F555W 


u2j20d 


1380 


F814W 


u2j20d 


830 


6.4 


NGC4649 


E2 


15.6 


2.5 






F555W 


u2qo03 


2100 


F814W 


u2qo03 


2500 


87. 


NGC5018 


E3 


42.6 


3.0 






F555W 


u3m725 


800 


F814W 


u3m725 


1800 


4.3 


NGC5845 


E 


24.0 


0.27 






F555W 


u30709 


2140 


F814W 


u30709 


1120 


3.8 


NGC5846 


EO-1 


21.1 


1.5 






F555W 


u36j04 


1300 


F814W 


u36j04 


1400 


5.4 


Note. — Summary of the sample galaxies and the observations used in the present analysis. The morphological type (Type) 



is taken from NED. All distances ^Disij axe ^ 



the Cepheid zero-point (e.g. Jensen et al. 2003), except that of NC 



extrapolated K s -band luminosities [Lk ) were taken from 2MjAS2-^ 
chosen for consistency with the stellar popi 



f M 



araston 



corrected by -0.16 mag 
wVinr-Vi was taken from 



tr> account 



visions to 
). Total, 



"For - each 



magnitudes do not differ significantly (e.g. Carpenter 
(Obsld) of the WFPC2 datasets used 
the approximate semi-major axis of the central, excluded region (R m in). 



Fabcr et al. (h! 
opting Mkq — - 3.41. 1 he zero-point was 
) and assuming the K-band and K s band 
>V filter band, we list the observation ID 
the actual filter used (Filt) and the total exposure in seconds (Exp). In addition, wc show 



partly overlaps the image of NGC 4649, and so to mini- 
mize confusion wc excluded any sources lying within the 
B-band twenty-fifth magnitud e ellip se of NGC 4647, as 



listed in de Vaucouleurs et al. (1991) 



Photometry was obtained with a 3 pixel radius aper- 
ture for the PC, and a 2 pixel aperture for the WF CCDs. 
The background was accumulated from an annulus be- 
tween 3 and 5 pixels in radius. Since a GC may be 
slightly resolved in the case of the galaxies here, the prox- 
imity of the background annulus to the centre of the GC 
may lead to a slight over-subtraction and hence underes- 
timate of the true brightness of a large GC. To estimate 
the likely magnitude of this ef fect, we sim ulated images 
of a GC at 15 Mpc, adopting a King| (1962) model for its 
surface brightness profile. We adopted a "concentration" 
(the logarithm of the ratio between the truncation and 
core radii) of 1.25 and a half-light radius of 2.5 pc (e.g. 



ti qn correction, we con verted the line-of-sight estimates 
of ISchlcgcl et al.| (|l998| , listed in NED) to the wavelength 
of the mid-poi nt of the appropriate filter, using the ex- 
tinction law of I Cardelli et alj ( 1989] ). This was typically 

u. value 1 



very close to the Schlcgcl et ai. value for the nearest stan- 
dard photometric band. The photometry was converted 
to the standard U BVPJ system foll o wing t he conversion 



formulae given in Holtzman et al. (1995a), assuming a 



standard V-I colour of 1.14 and a B-V colour of 0.97 
For a consistent comparison between the galaxies, we 
converted the (B,I) band photometry of NGC 1399 and 
NGC 1404 into (V,I) band p hotometry by adopting the 
conver sion formula given in |Gcbhardt fc Kisslcr-Patig 
Q1999D. 



[Kundu et al 



1999). The images were convolved with im 
ages of the PSF for the wide field and planetary cameras, 
generated with the TinyTim 3 software package. For the 
planetary and wide field cameras, the magnitude of over- 
subtraction was ~0.1 mag and ~0.02 mag respectively. 
This magnitude is similar to the <0.1 ma g error on the 



t otal magnitude of the GC estimated by Kundu et al 



(1999) on account of typical GCs being slightly resolved 
at ^ 10-20 Mpc. These small effects will not significantly 
affect the interpretation of our results. Photomet ric zero- 
points were taken from Holtzman et al. (1995a), taking 
into account the appropriate gain correction. Aperture 
correction factors were estimated based on simulations 
using TinyTim for each CCD, assuming the spectrum 
of a K4V star, and generally agreed with the calibra- 



tion of Holtzman et al. ( 1995b| ). To estimate the extinc- 



3 Available at http://www.stsci.edu/software/tinytim/ 



3. GC COLOURS 

3.1. Colour Histograms 

In order to investigate the colour distribution of GCs, 
we computed histograms of V-I colour for each galaxy. In 
Paper I we investigate the dependence of the probability 
that a GC hosts an LMXB as a function of metallicity, 
which we estimate by a simple transformation of the GC 
colour. Although, in general the colour distributions of 
GCs in early-type galaxies are bimodal, we do not expect 
whether a GC belongs to the red or blue subpopulation 
to affect its likelihood to host an LMXB, but rather this 
should depend on the metallicity. To obtain a global fit 
to this metallicity dependence in Paper I, we therefore 
need only to parameterize the composite colour distri- 
bution of all the GCs with a model which describes its 
overall shape reasonably well, and we do not need to cap- 
ture the fine detail. As a simple approximation, we chose 
to adopt a single Gaussian m odel, which has previously 
been fitted to such data (e.g. Gebhardt fc Kissler-Patig 
1999). We show in Fig [I] the Gaussian fit to the com- 



4 



Humphrey 



1 1.2 

V-I (mag) 



Fig. 1. — Composite GC colour distribution, fitted with a sin- 
gle Gaussian model. Although the model does not capture all of 
the fine detail of the distribution (which may be bimodal), it is a 
reasonable approximation of the overall shape. 

posite GC colour histogram for all of the galaxies, with 
a mean V-I of 1.033±0.005 and a of 0.186±0.005. Al- 
though the fit is formally unacceptable (x 2 /dof=57/27), 
as is clear from Fig [j] the model does describe the global 
shape of the distribution reasonably well; in fact it is 
accurate to ~15% for most of the interesting range of 
colours. This is adequate for our purposes in Paper I, as 
evinced by the excellent agreement we find in that paper 
between our measured metallicity dependence and that 
found by other authors with different techniques. 

For comparison with other authors' work, and to fa- 
cilitate the comparison between galaxies in Paper I, we 
also fitted the single Gaussian model to the GC colour 
distributions of each galaxy. Given the fewer numbers 
of sources in each bin, we adopted a Cash-C minimiza- 
tion algorithm to fit the model to the data, and assessed 
the goodness-of-fit via 100 Monte-Carlo simulations. For 
each simulation, an artificial dataset was created given 
the best-fitting number of sources in each colour bin, 
which was then fitted with our preferred model. The 
fraction of simulations yielding a higher Cash-C statis- 
tic (i.e. worse fit) than the model fit to the real data 
can be used as an estimate of the goodness of fit. The 
fit was formally acceptable at a 5% significance level in 
most of the galaxies. Unacceptable fits were found for 
NGC1399, NGC3115, NGC4494, NGC4552, NGC4621 
and NGC 4649, which typically reflects evidence of bi- 
modality or asymmetry in the GC colour distribution. 
Nonetheless, the fits of the single Gaussian model to the 
GC colour distributions of each galaxy (shown i n Table [ 



gene rally yield r esults consistent with those of Gcbhardt 
fc Kikslcr-Patip (119991 ), [Kundu fc Whitmorcj pX)Ta| Jb|) or 



theirs. 



Larsen et al.~[2001), where our sample overlapped with 



-3t2t — Mass - mctallicity relation 



We show in Fig |2j (left panel) the composite colour- 
magnitude diagram for all of the GCs detected in more 
than one filter. We also show representative error-bars 
and the data rebinned into a series of narrow magnitude 
bins. The data do not show convincing evidence of a sig- 
nificant colour-magnitude correlation. For example, the 
binned-up data are consistent (x 2 /dof=13.6/8) with no 



variation in colour with magnitude. The fit does not im- 
prove significantly (x 2 =12.2/7) when fitted with a linear 
model of the form {V-I) = (1.043±g;gJ|) + (0.003 ± 
0.005)V; this is consistent with no variation within er- 
rors and indicates that, if any correlation is present, it is 
extremely weak. The V-I colour used in the present work 
does not cleanly separate the metal-poor and metal-rich 
sub-populations of GCs; coupled with the small field of 
view of the WFPC2, which limits the numbers of blue 
GCs detected, we do not see the so-called "blue tilt" (i.e. 
colour-magnitude correlation) recently observed in the 
metal-poor GC populations o f several galaxies (Harris 



eTaT]|200C[: |Strader et aLlpOOj ) 

In the right panel ot Fig 



we show the same data 
converted into mass and metallicity. To achieve this, 
we linearly interpolated the V and I-band photometry 
onto the updated s imple stellar p opulation (SSP) model 
grids of Maraston ( 2005 , see also Maraston 1998 ) 4 . We 
adopted these models for consistency with our previous 
study of the mass and metallicity of early- type galaxies, 
in wh ich we adopted the al pha-enhanced Lick index mod- 
els of Thomas et al 
(Humphrey & Buotc 
profiles of early-type galaxies 
we found that our measured 



(120031) 



built upon these SSP models 
2006|) . I n our study of the mas s 

" fg§ ), 



Humphrey et al 
band mass-to 



K 



igEF 



ra- 
tios in the (baryon-dominated) centres of th e galaxies 
were consistent with the predictions of t he [Maraston 
models, provided one adopted a Kroupaj (2001) initial 
mass fu nction, wh ich we therefore adopted in the present 
work. Maraston provides the mass-to-light ratio as a 
function of metallicity and age. We assumed an age 
of 13 Gvr and adopted the SSP models computed for 
a Kroupa initial mass function and a blue horizontal 
branch (HB) for the GC stars. A small fraction of the 
GCs with best-fitting abundances pegged at the high or 
low boundaries ([Z/H]=-2.25 and 0.67, respectively) were 
discounted. Since the V and I-band mass-to-light ratios 
show a dependence upon the metallicity, the lack of a 
correlation between V-I and V actually translates into 
a slight mass-metallicity correlation. Formal correlation 
testing was carried out using three different methods, 
Pearson's linear correlation test, Spearman's rank-order 
test and Kendal's r test. All three tests indicated that 
the probability of no correlation is less than 10 -29 . 

We also show in Fig || the data binned into a 
series of narrow mass bands. The error-bars on 
the mean GC metallicity in each bin were estimated 
from the 1-a scatter of the data-points, using the 
central limit theorem. These data are well-fitted 
(x 2 /dof=4.3/5) by a simple linear model of the form 
[Z/H] = (-2.1 ± 0.2) + (0.25 ± 0.04)logi M, where M is 
the mean GC mass. To estimate the intrinsic scatter 
about this best-fitting relation, we additionally fitted this 
model to the y-on-x data-points for each GC, taking into 
account the errors on both axes with a n algorit hm similar 
to that out lined in Press et al. ( 1992 , see also Trcmainc 
2002] ) . We added in quadrature with the metallic- 



et al 



ity error an additional, constant, systematic error-term, 
which was adjusted until the x 2 /dof value was 1, at which 
point its value is approximately equal to the 1-a intrinsic 
scatter about the best-fitting linear relation. In this way, 

4 http: / /www-astro. physics. ox. ac.uk/~maraston/Claudia's_Stellar 
_Population_Models.html 
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v-i 



log 10 (M/M o) 



Fig. 2. — Left: Composite colour-magnitude diagram for all of the GC candidates. We show (blue) representative error-bars for a range 
of V magnitudes. Overlaid we show the data-points binned into a series of magnitude bins, clearly indicating little or no colour-magnitude 
relation. Right: the same data, transformed into mass and metallicity space (see text). Overlaid we show the same data binned into a series 

of metallicity bins, clearly indicating a slight mass-metalltcity r-nrrplatir.Ti Ji^ also show (dotted line) an extrapolation of the best-fitting 

mass-metallicity relation for early- type galaxies found by Thomas et al. (200. 



X 




element mass fraction ("Z") and [Z/H] as those authors. 
The resulting mass-metallicity plot is shown in Fig 0. We 
also show the data binned into a series of narrow mass 
bands, which clearly indicates a mass-metallicity correla- 
tion similar to that found above. Fitting the binned data 
with a simple linear model, we found that the data are 
inconsistent with a constant [Z/H] model (x 2 /dof=42/6), 
but are in agreement (x 2 /dof=ll/5) with a model of the 
form [Z/H] = (-1.8 ± 0.2) + (0.12 ± 0.03)logi M. Both 
the slope and normalization of the bes t-fitting mo del dif- 
fer from the relation obtained with the Maraston models, 



l O g 10 (M/M o) 

Fig. 3.— Same as Fig|, right panel, but using the PEGASE SSP 
models to convert colour and magnitude to mass and metallicity. 

we estimated this scatter to be 0.62 dex in metallicity. 

Since the assumption of a blue HB is most appropri- 
ate for metal-poor GCs, we have investigated the impact 
of using the red HB models instead of the blue models. 
We found that this produces a small shift in [Z/H], at 
maximum ^0.2 dex, which does not appear to correlate 
with mass. We found similar shifts (albeit typically of 
a different sign) if we adopted a GC age of 10 Gyr. We 
therefore conclude that relaxing our assumptions of a 
fixed GC age or a blue HB for all GCs is unlikely to 
make the observed trend disappear. In order to ver- 
ify that our results do not sensitivity depend upon the 
choice of SSP models, we also experimented with adopt- 
ing models based o n the PEG ASE. 2 code of Fioc & 



Rocc a- Volmerange (1997, 199E) to perform the transfor- 
mation from colour to mass and metallicity space. We 
t ook most of the default parameters of the code, adopted 



a Kroupa IMF, assumed a single burst of star- formation 
and wc included the mass of neutron-stars and white 
dwarfs in the calcula tion of the mass-to-light ratio. For 



consistency with the Maraston models, we adopted ap- 
proximately the same transformation between the heavy 



which reflects quantitative differences in the relation be- 
tween the mass-to-light ratios and the metallicity for the 
two sets of models. 

4. LUMINOSITY FUNCTIONS 

4.1. Source Detection Incompleteness 

In order to measure accurately the GCLF in each 
galaxy it is essential to take into account the source 
detection incompleteness. The ability of the algorithm 
outlined in § § to detect point-sources is a strong func- 
tion both of the source flux and the "background" level, 
which varies dramatically over the WFPC2 field of view. 
It is therefore certain that some fraction of the GCs in 
any given galaxy will not be detected. To estimate the 
completeness of our source lists, we therefore carried out 
extensive Monte-Carlo simulations. Starting with the 
V-magnitude (B-band for NGC 1404 and NGC 1399) of 
the faintest sources detected in each galaxy, we stepped 
through a series of magnitude bins, each separated by 
0.2 mag, in each of which we performed 10 Monte-Carlo 
simulations to estimate the completeness level. We have 
explicitly confirmed that starting from one magnitude 
fainter than this limit makes no measurable difference to 
our key results for one galaxy, NGC 4649. Each simula- 
tion entailed the addition of ten simulated point sources 
to the actual WFPC2 images of the galaxy in each filter. 
The V-I colours of the sources were chosen randomly by 
drawing, with replacement, from the V-I distribution of 
the detected sources. The source detection algorithm was 
run on each image and the number of simulated sources 
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TABLE 2 

Derived Properties of the GC Populations 



Name 


]\rdet 
JV GC 




V-I 




Ov-I 


f. GC 


L GC 


r FOV 

(10 lo L G 




T 

m v 






S N 




Sl 
















(10 7 L Q ) 


















NGC1332 


198 





.998 ± 0.03 


0. 


un D3 

z -0.02 


0. 


.60 


8.3 ± 1.2 


2.0 ±0. 


9 


24. 


,4±0. 


A 


2 


.5 ± 1.1 





.42 ± 0.19 


NGC1387 


28 











.22 


0.79 ±0.51 


1.2 ±0. 


3 








0. 


.38 ±0.26 


0. 


.065 ± 0.045 


NGC1399 


467 


1 


.034 ± 0.009 


0. 


OQ +0.007 
1^3_ 0Q6 





.81 


12. ± 1. 


2.0 ± 1. 


.8 


24 


,0±0 


.2 


3 


.5 ± 3.1 





.59 ± 0.53 


NGC1404 


109 


1 


.01 ± 0.02 


0. 


.14 ± 0.01 





.74 


5.3 ± 0.9 


1.7 ± 1. 


.8 


24 


,3±0 


.1 


1 


.8 ± 2.0 





.32 ± 0.35 


NGC1553 


71 





.996 ± 0.04 


0. 


i ~+0 04 
la -0.03 





.76 


2.1 ± 0.5 


2.7 ±0. 


6 


24 


,5±0 


.6 





.45 ± 0.14 





.077 ± 0.024 


NGC3115 


136 


1 


.01 ± 0.02 


0. 


1 £+0.02 

15 -o!oi 





.92 


2.3 ± 0.4 


1.1 ± 0. 


1 


22 


,7±0 


.3 


1 


.2 ± 0.3 





.21 ± 0.04 


NGC3585 


90 





.99 ± 0.04 


0. 


QQ+0.05 

^-0.03 





.83 


2.4 ± 0.4 


3.0 ± 0. 


.8 


24 


6±0 


.4 





.47 ±0.15 





.081 ± 0.026 


■\T /~i ocn7 

1N(jC3dU7 


106 





.95 ± 0.03 


0. 


1Q +0.03 
1M -0.02 





.62 


4.9 ± 0.9 


2.4 ± 0. 


.8 


24 


.0±Q 


.4 


1 


.2 ± 0.4 
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.20 ± 0.07 


NGC4125 


196 


1 


.05 ±0.03 


0. 


24 +0.04 
z -0.03 





.74 


2.6 ±0.4 


3.2 ±0. 


.8 


25 


,5±0 


.3 


1 


.3 ±0.5 





.082 ± 0.023 


NGC4261 


241 


1 


.00 ±0.02 


0. 


.20 ±0.02 





.56 


11. ± 1. 


3.4 ±3. 


.4 


24 


,6±0. 


.3 


2 


0±2.0 





.34 ± 0.34 


NGC4365 


291 


1 


.01 ±0.02 


0. 


.18 ±0.01 





.81 


8.2 ±0.9 


2.3 ±0. 


6 


24 


,0±0 


.2 


2 


.1 ±0.6 





.36 ±0.11 


NGC4472 


397 


1 


.08 ± 0.02 


0. 


.18 ±0.01 





.83 


7.1 ± 0.7 


3.9 ±0. 


5 


23 


,8±0. 


.2 


1 


.1 ±0.2 





.18 ±0.03 


NGC4494 


155 





.98 ± 0.02 


0. 


.18 ±0.02 





.86 


3.0 ±0.4 


1.5 ±0. 


5 


23 


,5±0. 


.2 


1 


.2 ±0.4 





.20 ±0.07 


NGC4552 


240 


1 


.02 ±0.02 


0. 


-ia+0.02 





.85 


4.2 ± 0.5 


1.3 ±0. 


1 


23. 


4±0. 


.2 


1 


.9 ±0.6 





.32 ±0.10 


NGC4621 


173 


1 


.04 ±0.02 


0. 


ic+0.02 
lo -0.01 





.86 


3.9 ± 0.5 


1.7 ±0. 


5 


23. 


4±0. 


.2 


1 


.3 ±0.4 





.23 ±0.07 


NGC4649 


323 


1 


.05 ±0.02 





.17 ±0.01 





.87 


6.8 ± 0.7 


2.9 ±0. 


1 


23. 


,6±0. 


2 


1 


.4 ±0.3 





.23 ±0.04 


NGC5018 


90 





S '-0.05 


0. 


01 +0.05 
zl -0.03 


0. 


.47 


7.8 ± 1.8 


9.7 ±4. 


1 


25. 


1+0.6 

4 -o.5 







.46 ±0.22 





.080 ± 0.038 


NGC5845 


45 





MD -0.15 


0. 


nn+0.19 
ZM -0.07 


0. 


.47 


1.4 ±0.8 


1.9 ±1. 


1 








0. 


.45 ± 0.36 


0. 


.077 ±0.062 


NGC5846 


288 


1 


.04 ±0.02 


0. 


-ia+0.02 





.76 


5.3 ±0.6 


2.6 ±1. 


.1 


24. 


,8±0 


.2 


1 


.4 ±0.6 





.21 ±0.09 



Note. — Derived properties of each of the galaxies' GC distributions. We list the number of GCs detected (Nf^), the 
mean colour of the GCs (V-±o), the half-width of the V-I distribution (cry_j), the total fraction of the GC luminosity 
which has been detected (/^ ), the total V-band luminosity of the GCs, correcting for incompleteness {Lgc)-, the 
V-band magnitude of the galaxy in the field of view (L^ OV ), the V-band GCLF turnover (m^), the specific frequency 
of GCs (Sn) and the specific luminosity (Sl)- The correlations between these properties and the LMXB distribution are 
discussed in detail in Paper I. All quoted errors arc-90% confidence limits, incorporating only statistical uncertainties. 
An assessment of statistical errors is given in Table pi. 



TABLE 3 

Systematic error-budget for derived parameters 



Name 


fGC 




Lgc 












Sl 
















A/ ree 


^■wid 




^free 


^wid 


A s t a t 


^free 






NGC1332 


-0.052 


-0.063 


1.2 


0.83 


0.98 


1.1 


1.5 


1.4 


0.19 


0.042 




0.050 


NGC1387 




-4.8 x 10" 4 


0.51 




3.9 x 10" 3 


0.26 




0.71 


0.045 






3.2 x 10" 4 


NGC1399 


-0.021 


-0.049 


1.0 


-0.50 


0.88 


3.1 


0.97 


1.0 


0.53 


-0.025 




0.043 


NGC1404 


-0.068 


-0.055 


0.86 


-0.18 


0.36 


2.0 


0.77 


1.3 


0.35 


-0.011 




0.021 


NGC1553 


-0.18 


-0.18 


0.48 


0.14 


0.20 


0.14 


0.74 


0.91 


0.024 


5.0 x 10 


-3 


7.4 x 10" 3 


NGC3115 


-0.022 


-0.039 


0.36 


-0.46 


-0.14 


0.23 


0.29 


0.32 


0.040 


-0.040 




-0.012 


NGC3585 


-0.095 


-0.11 


0.44 


-0.42 


-0.23 


0.15 


0.31 


0.47 


0.026 


-0.014 




-7.7 x 10- 3 


NGC3607 


-0.016 


-7.6 x 10~ 4 


0.93 


0.29 


0.61 


0.43 


0.21 


0.18 


0.074 


0.012 




0.025 


NGC4125 




-0.089 


0.37 




0.27 


0.53 




0.61 


0.023 






8.5 x 10" 3 


NGC4261 


-0.11 


-0.032 


1.5 


3.7 


1.7 


2.0 


2.6 


-0.023 


0.34 


0.11 




0.050 


NGC4365 


-0.043 


-0.016 


0.86 


0.044 


0.80 


0.61 


0.59 


0.59 


0.11 


1.9 x 10" 


-3 


0.035 


NGC4472 


-0.035 


-0.042 


0.69 


-1.3 


0.051 


0.16 


0.55 


0.30 


0.028 


-0.034 




1.3 x 10" 3 


NGC4494 


-3.5 x 10~ 3 


-0.085 


0.42 


-0.32 


0.95 


0.42 


0.41 


0.053 


0.071 


-0.022 




0.065 


NGC4552 


-0.032 


-0.060 


0.51 


-0.61 


0.66 


0.57 


0.95 


0.20 


0.098 


-0.046 




0.051 


NGC4621 


0.011 


-0.014 


0.51 


0.87 


1.8 


0.40 


-0.039 


-0.066 


0.069 


0.051 




0.11 


NGC4649 


-0.017 


-0.014 


0.68 


0.22 


1.5 


0.25 


0.19 


0.23 


0.041 


7.6 x 10 


-3 


0.050 


NGC5018 


-0.011 


0.097 


1.8 


1.2 


0.043 


0.22 


0.10 


-0.13 


0.038 


0.013 




4.4 x 10" 4 


NGC5845 




-0.019 


0.59 




0.59 


0.43 




7.4 x 10" 3 


0.073 






0.045 


NGC5846 


-0.014 


-0.048 


0.57 


-0.74 


-0.090 


0.62 


0.64 


0.61 


0.088 


-0.029 




-3.5 x 10~ 3 



Note. — Error-budget for fj , Lgc, Sn and Sl, showing the statistical errors (A. s tat) and an estimate of the impact on the 
parameters from allowing the GCLF turnover to fit freely while fitting all of the available data (Af ree ) or from adopting a GCLF width 
of 1.5 mag (A. w id)'-, see text. We indicate by ellipses cases where the data could not be constrained (NGC 1387 and NGC 5845) or where 
the parameters were already derived allowing the GCLF turnover to be free (NGC 4125). 
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FlG. 4. — Comparison of the measured GCLF n turnover and 
the corrected I-band SBF distance moduli (Table |l|). NG Q5Q1£ 



for which the distance was derived from the D n -<r relation (Faber 
et al.l 198E), is shown as a diamond. The straight line denotes 
the best-fitting straight line relation (see text), which is consistent 
with previously-reported correlations. This demonstrates the util- 
ity of the GCLF turnover magnitude as an approximate "standard 
candle" . 

which were detected, as well as their measured magni- 
tudes, was determined. Once 98 percent of the simu- 
lated sources were detected in any bin, the simulations 
were stopped and all sources brighter than the corre- 
sponding magnitude were assumed detected at their true 
magnitude (this is, in general, a very good approximation 
since the measurement errors on the brightest GCs are 
typically much less than the 0.2 mag binsize we adopt). 
The sources were assumed to be distributed as the opti- 
cal light, which, for these purposes was modelled as an 
elliptical de Vaucouleurs profile centred at the peak of 
th e galaxy emission, with the size and orientation given 
in d e Vaucouleurs et aL ( 1991, except for NGC 1332, for 



whic i the position angle was modified to match better 
the observed image of the galaxy). Point-source images 
were generated with the TinyTim tool for the appropri- 
ate filter and focal plane position of the source, assuming 
the spectrum of a K4V star. 

Any given parameterization of the GCLF which we 
fit predicts the number of sources which should, in the 
idealized case, be observed in any given magnitude bin. 
To account for incompleteness and statistical errors, this 
model was convolved, in each bin, with a function, de- 
termined from our simulations, which redistributes the 
sources into adjacent bins, based on the likelihood it is 
detected in each magnitude bin. The corrected model 
can then be fitted directly to the data. This method 
is ex actly analogous to computing a "response matrix" , 



of magnitude bins of width 0.2 mag. We created an asso- 
ciated response matrix, as described in Appendix [A] and 
fitted the data with dedicated software through a Cash-C 
statistic function minimization procedure. We modelled 
the GC magnitude distribution as a simple Gaussian, i.e. 

dn GC N TOT ( (V-V T ) 2 \ n , 

^ = vW ex H (1) 

where ncc is the number density of GCs, as a function of 
V-band luminosity (V) , Ntot is the total number of GCs 
and Vt is the apparent V-band peak luminosity ( "GCLF 
turnover"). This has been previously shown to be a rea- 



sonable description of the G CLF shape (e.g. Harris 1991 ; 
Kundu & Whitmorc 2001a). For consistency with past 
analysis, and to minimize the sensitivity of our results 
to the our source detection incompleteness algorithm, 
we initially fitted only those data-bins for which the 
source detection completeness was estimated to be 50% 
or greater. In general, our incompleteness cut allowed 
insufficient data below the GCLF turnover to constrain 
all three para meters of the Gaussian fit a nd so, based 



on pa st work (Kundu & Whitmorc 2001a; Jordan et al 
2006[) , we fixed the GCLF hall- width (cry) to 1.3 mag 
This should be a reasonable average value for giant ellip- 
tical galaxies, although in individual objects significantly 
broader or narrower GCLFs may exist. 

In Table ^|we list the V-band turnover magnitude, my 
for each galaxy. To enable a consistent comparison, we 
convert B-band turnover magnitudes for NGC 1399 and 
NGC 1404 to the V-band by assuming a constant B-V of 
0.97, consistent with the colours of these galaxies, and 
also consistent with a ~ solar metallicity 10 Gyr stellar 



population ( Maraston| |2005| ) . We ignored NGC 1387 and 
NGC 5845, lor which there were too few measured GCs 
to make the GCLF turnover measured in this way re- 
liable. Where our samples overlap, our mea sured my 



were in go o d agre ement with those fo u nd by Kundu 
Whitmorc! fl2001a| Jfcl) and |Larscn et aT| (|200l| , who fitted 
a slightly different functional form to the GCLF), giving 
us confidence in our fitting procedure. 

In Fig | we show my as a function of the distance 
modulus. Taking into account errors on both the x 
and y data values we fitted the data with a function of 
the form my = (m — M)sbf + My, where My, the ab- 
solute turnover magnitude, was fitted freely. The fit 
was formally unacceptable (x 2 /dof=35/16), indicating 
significantly more scatter about this relation than the 
statistical errors. Omitting the most discrepant data- 
point, that of NGC 4125, significantly improved the fit 
(x 2 /dof=23/15), although it remained only marginally 
acceptable. For this fit we constrained M^=— 7.3 0±0.05, 
in good agreement with previous estimates (e.g. Kundu 



thro ugh which any parameterized model can be folded in fc Whitmore] |2001a|) , especially when we ac count for re- 



X-ray spectral-fitting (e.g. Davis " 2001 ). We discuss this 



method, and the intrinsic uncertainty associated with our 
method of computing the matrix in Appendix |A| This 
procedure not only accounts for source detection incom- 
pleteness, but also the Eddington bias, while preserving 
the statistical integrity of the data. 

4.2. Fitting the GCLFs 

For each galaxy, we accumulated a V-band (B-band 
for NGC 1399 and NGC 1404) GCLF by constructing a 
histogram of the numbers of sources observed in a series 



cent r evisions to the Cepheid zero-point (e.g. Jensen et al 
2003). In general, however, the ad hoc removal of data- 
points should be avoided. Taking our fit to NGC 4125 
at face valu e im plies intrinsic scatter (estimated as de- 
scribed in § 3.2) of 0.4 mag, with a 90% lower limit of 
0.26 mag, and a best-fitting of -7.26 ± 0.15. 

In order to compare the LMXB and GC populations 
in Paper I, we require an estimate of the fraction of the 
GC luminosity which is detected in the WFPC2 field of 
view, f- ; . To estimate this for each galaxy, we adopted 
the canonical GCLF (My=-7.3, ov=1.3) to each dataset 
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(except for NGC4125, for which we adopted My de- 
rived from fitting) and computed the observed and to- 
tal luminosity of GCs by appropriately integratingthc 
incompleteness-corrected and the uncorrected models. 
We tabulate ff c in Table |. 

4.3. Specific frequency and luminosity 

To quantify the richness of the GC populations of each 
galaxy within the WFPC2 aperture, we computed two 
quantities, the GC specific frequency, S n. and the GC 
specific luminosity, Sl (e.g. Harris 1991 ), defined by the 
relations: 



5' 



JV 



7V GC 10°- 4 ( M " OV+15 > 



S L = 100L GC /L 



FOV 
V 



(2) 
(3) 



where Nqc is the total number of GCs inferred from fit- 
ting the GCLF, Lqc is the total luminosity of the GCs 
inferred from the fit, My ov is the V-band absolute mag- 
nitude of the galaxy in the field of view and Ly OV is 
the corresponding luminosity. Although we find signifi- 
cant intrinsic scatter in My, in order to produce inter- 
esting constraints, we computed Nqc and Lqc having 
fixed My to -7.3 for each galaxy, except for NGC4125, 
for which we allowed the turnover to fit freely since it 
is substantially fainter than My (Fig f4j). We estimated 



M^ ov , following |Kundu fc Whitmorc| (p001a|), by inte 



grating the total counts in the field of view (having ex- 
cluded obvious bright interloper sources) and adopting 
a background level of 0.01 and 0.052 e~ s~ 1 pixel~ 1 for 
the PC and WF CCDs. We computed these quantities 
in the V-band for all galaxies, excepting NGC 1399 and 
NGC 1404, for which they were computed in the B-band. 
For these two systems, we subsequently corrected Sn to 
the V-band, assuming a galaxy B-V colour of 0.97. The 
best-fitting results are shown in Table |[ Where our data 
overlap, our computed Sn (shown in T able | |) generally 
agree very w ell with those reported by |Kundu fc Whit- 
" fl2001a| ,|b|). The GCLF data and the best-fitting 



lown in Fig |^. 



moclcls are st 

This model fits the data brighter than the cut-off well. 
However, extrapolation to fainter magnitudes generally 
reveals some discrepancies (Fig fsj). This may reflect ac- 
tual differences between the GCLF width and our canon- 
ical ay = 1.3 mag, the increasing importance of (and diffi- 
culty in eliminating) interlopers at faint magnitudes, and 
possible deficiencies in our incompleteness-estimation al- 
gorithm. To investigate how important these effects may 
be, we experimented with fitting all of the data-points 
and freeing my. The best-fitting models for such fits 
are shown (dashed lines) in Fig |[ In Table ||, we tab- 
ulate the systematic uncertainties on Sn and Sl arising 
from this choi ce o r, similarly, fixing the GCLF width to 
ov=1.5 (see § 4.4). Another potential source of system- 
atic uncertainty is the assumed spatial distribution of the 
point sources in our simulations. In order to assess how 
significantly this may affect our results, we recomputed 
the incompleteness correction for one system, NGC 4649, 
but distributing the simulated sources as a (3 model, with 
[3 = 0.46 and core radius 0.81'. Such a model approxi- 
mately fits the distribution of the LMXBs in this system. 
We found Sn changed by 0.20 and Sl by 0.026, which are 
consistent with the statistical errors, and other sources of 
systematic uncertainty for this object. In general, since 



Sl is more sensitive to the bright end of the GCLF dis- 
tribution, we found it was more robust than Sn to these 
choices; the systematic effect on Sn of freeing My was 
generally larger than the statistical errors, whereas the 
impact on Sl, while not negligible, was typically smaller 
than, or comparable to, the statistical errors. Sl also has 
the advantage that it is independent of the any error in 
the assumed distance to the galaxy. 

4.4. Composite GCLF 

In Paper I we investigate the dependence of LMXBs in 
GCs by comparing the composite luminosity function of 
the GCs in all galaxies with that of the GCs containing 
LMXBs. In order to investigate the composite GCLF, we 
first computed V-band GCLFs for each galaxy in a set 
of standard absolute magnit ude b ins, and appropriate re- 
sponse files, as described in § 4.1. To mitigate against the 
intrinsic scatter in My, for each galaxy we adopted a dis- 
tance modulus given by the measured my+7.3 to convert 
from apparent to absolute magnitude. To take account 
of incompleteness effects the response matrices were av- 
eraged using the Heasoft tool addrmf weighting each ma- 
trix by the number of detected sources. To enable us to 
combine the GCLFs for NGC 1399 and NGC 1404 with 
the other data, we converted from the B to V band by 
assuming a constant B-V=0.97. 

The combined GCLF was well-fitted by a single Gaus- 
sian model, with cry = 1.51 ± 0.04. The data and best- 
fitting model are shown in Fig |[ We note that the best- 
fitting Gy is slightly larger than our adopted 1.3 mag 
used in computing Sn and Sl- This most likely arises 
due to the GCLFs of some, although not necessarily all, 
of the objects being slightly broader than we have as- 
sumed. The nature of a composite GCLF is that its 
shape is heavily biased towards the most GC-rich galax- 
ies, and so ay =1.5 need not be repr esentative of the ma- 
j ority o f the galaxies in our sample. Kundu & Whitmorc 
( |2001a|) found that NGC 4472 and NGC 4649, which ac- 
count for ^20% of the detected GCs, may have ay> 1.5. 
The GCLF-determined distance moduli have statistical 
uncertainties typically of the order ^0.1-0.2 (1-u), which 
should also cause the measured composite GCLF to be 
slightly broadened. Nonetheless, we estimate in Table || 
the systematic effect on the best-fitting Sn and Sl of 
fixing cry=1.5. 

5. SPATIAL DISTRIBUTION OF GLOBULAR CLUSTERS 

In order to investigate the spatial distribution of the 
GCs, we separately binned radial distributions for the 
red and blue GCs, defining any cluster with VT>1.1 as 
red, and the remaining clusters as blue. We note that 
this does not cleanly separate the blue and red sub- 
populations of clusters and so measured differences in 
the spatial distributions probably underestimate intrin- 
sic differences between the red and blue sub-populations. 
However, it does divide clusters on the basis of their ac- 
tual metallicity (rather than asking to which subpopula- 
tion they belong), thereby allowing us to assess whether, 
for our galaxies, more metal-rich clusters are, on average, 
distributed differently to more metal-poor ones. This dis- 
tinction is extremely important when comparing to the 
spatial distribution of LMXBs (Paper I). 

To compute the radial profile, the radial distances of 
each source from the optical centre of the galaxy, identi- 





Fig. 5. — GCLFs of each of the galaxies in the sample, shown with the best-fit Gaussian model (solid line), corrected for source detection 
incompleteness. Data-bins for which source detection completeness was below 50% (not included in the fit) are drawn with dotted lines. 
The dashed line shows the best-fitting model if all of the data are fitted and the turnover magnitude is allowed to fit freely. Although the 
data were fitted by minimizing the Cash-C statistic, we show error-bars estimated from Gchrcls' approximation (for the upper limit) to 
guide the eye. In most cases the data are well-represented by the simple Gaussian model, and the results are not highly sensitive to the 
inclusion or exclusion of the more sources at the uncertain faint end of the GCLF. 



fied by eye, were binned into a set of 7 radial bins, with 
outer annuli 0.33, 0.66, 0.82, 1.0, 1.2, 1.3 and 1.5', re- 
spectively. The data were rebinned to ensure at least 3 
sources per bin and the radial distribution was fitted with 
a variety of models by a Cash-C statistic minimization 
algorithm. To correct the spatial distributions of sources 
for detection incompleteness, we computed a GCLF cor- 
rection function, similar to that described in § 4.1 for 



each bin. Assuming the GCLF does not vary spatially, 
this correction was combined with the measured GCLF 



for that galaxy to estimate the fraction of sources which 
would be detected in each bin. 

We show the spatial distributions of a sub-sample of 
the galaxies, with sufficient measured GCs to yield in- 
teresting constraints, in Fig @. We compare in this fig- 
ure these distributions to the optical light, as parameter- 
ized by a de Vaucouleurs model, with its effective radius 
fixed to that measured in the K-band f rom the 2MASS 
Extended Source Catalog ( Jarrett 200C ) . We found that 
both the red and blue GC distributions were significantly 
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Fig. 6. — Composite GCLF for all the galaxies in the sample. 
Also shown is the best-fitting single Gaussian model, corrected 
for source detection incompleteness. Although the data are fitted 
with the Cash-C statistic, to guide the eye we show error-bars 
estimated from the Gehrels' approximation. This demonstrates 
the good quality of a Gaussian GCLF fit to the data. 



TABLE 4 
Radial distribution fits 



Galaxy 



a Blue 



NGC1332 
NGC1399 
NGC3115 
NGC3585 
NGC3607 
NGC4125 
NGC4261 
NGC4365 
NGC4472 
NGC4494 
NGC4552 
NGC4621 
NGC4649 
NGC5846 



-+0.2 
-0.1 

-0.73 ±0.23 

„+0.44 
-0.36 
,+0.73 
-0.52 
1+0.6 
-0.4 
.,+0.41 
-0.35 
-,+0.30 
-0.25 

-0.79 ±0.29 
-0.61 ±0.24 
-0.97 ±0.18 
-1.1 ± 0.2 

-o.ssi^M 

-0.42 ±0.28 
-1.0 ±0.2 



-1.5;! 

-0.73 

-0.54;! 

-0.7V 1 

-l.lj 

-0.26;! 

-0.9uj 



-0.98] 



a+0.21 
-0.18 

-0.59 ±0.17 

-0.088l°« 

-0.94iS;S* 

-1.1 ±0.2 

-0.81 ±0.23 

-0 Qi+ - 20 

-0 66+ - 20 

-0.21 ±0.27 
n n 9 +0.47 

-0.62 ±0.23 
-0.85 ±0.22 
-0.58 ±0.19 
-0.871°;- 



Note. — Results from a power-law 
fit to the radial distribution of the GCs 
within a radius of 1.5'. The power-law in- 
dex, a is given separately for the red and 
blue GCs. The distribution of blue GCs 
tends to be flatter than for red GCs. 

flatter than the optical light. These data could be fitted 
adequately by a single powerlaw model, 



dn , 



GC 



dr 



(4) 



where driac I dr is the number density of GCs as a func- 
tion of radius, r, and a is a parameter of the fit, as is the 
normalization N. The results of the fits are tabulated 
in Table [| and shown in Fig [| The GC distributions 
are considerably flatter than typicall y reported for stud 



ies with larger radial coverage (e.g. |Kisslcr-Patig| |1 997 



and references therein), but this simply reflects the flat- 
tening of the GC d istributions in the galaxy core (e.g. 



Bassino et al. 2006| ). In general, the red GCs are sig- 
nificantly more centrally-peaked than the blue GCs, al- 
though there were a few exceptions to this trend. This 
is probably a result of the size of the error-bars, the lack 



of full azimuthal coverage, the relatively small numbers 
of sources detected and, critically, the small radial range. 
In any case, on average, the index for the blue GC radial 
distribution is less negative by ^0.15. 

6. DISCUSSION 

In order to facilitate a direct comparison between 
LMXBs in early-type galaxies and possible GC candi- 
dates, we have presented a catalogue of GC candidates 
and key derived properties based on the HST WFPC2 
centrally-pointed observations of 19 nearby early-type 
galaxies. The complete source catalogue, including co- 
ordinates and photometry, is given in Appendix The 
total fraction of the GC light detected and GC specific 
luminosities are given for each galaxy in Table |[ In Pa- 
per I we combine these data with archival Chandra obser- 
vations of early-type galaxies to show that the majority 
of LMXBs must have formed in GCs. 

We self-consistently derived various properties of the 
GC populations, including fits to the GCLFs. We report 
specific frequency, Sn, and specific luminosity, Sl, of the 
GC populations. Since the luminosity is dominated by 
clusters brighter than the turnover, we found that Sl is 
much less sensitive than Sn to uncertainties in the GCLF 
turnover or width, making it a more robust indicator of 
the richness of the GC population where a significant 
fraction of the GCs are undetected due to source detec- 
tion incompleteness. We find that the absolute magni- 
tude of the GCLF turnover (My) exhibits intrinsic scat- 
ter from galaxy to galaxy of ~0.3-0.4 mag, when com- 
pared to SBF distance m o duli, co nsistent with the es- 
timate of Ferrarcsc et al. (2000). Kundu & Whitmorc 
( |2001a ) argued the intrinsic uncertainty in the GCLF 
distance measure is ~ 0.14 mag, but only when using 
a small subset of the galaxies in their sample with the 

However, there is a signif- 



smallest error-bars on my 
icant correlation between 



ly-nivmodci) and the error 
on my (as expected, since, the fainter the turnover, the 
less complete the data) , and so this systematically selects 
against galaxies with faint My, possibly reducing the 
scatter. Provided their error-estimates are reasonably 
r epre sentative, one can adopt the procedure outlined in 
§ |3.2| to estimate the intrinsic scatter for the 25 galax- 
ies m their full sam ple which overlap the SBF sample of 
Lbnry et al. (2001). Using this approach, we estimate an 
intrinsic scatter of ^0.38 mag (with a 90 % lower limit 



of 0.2 3 mag), c onsistent with our results. Jordan et al 
(200(3) (see also Jordan et al. 2007j ) fitted Gaussian mod- 



els to the GCLFs of early-type galaxies in Virgo, and 
found a trend of increasingly faint turnover magnitude 
with fainter Mb- For their galaxies with Mb<-20, we 
estimate by eye an intrinsic scatter of ~0.3 mag, also 
consistent with the estimate from our data. To some ex- 
tent, this scatter is driven by the presence of significant 
populations of diffuse star clusters in two systems; it may 
be that the presence of similar populations in the galax- 
ies we consider here have contributed to the significant 
scatter we observe. 

In keeping with past studies of GCs using V and I- 
band photometry we did not find any evidence of a sig- 
nificant correlation between colour and luminosity. Con- 
verting the colour and magnitude of each GC to mass and 
metallicity, however, we found evidence of a weak mass- 
metallicity correlation for the GC populations as a whole. 
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Fig. 7. — Azimuthally-averaged radial distribution of the GCs in a sub-sample of the galaxies. The radial distribution of the blue GCs 
are displaced upwards by an arbitrary amount for clarity. The measured radial profiles data-points (solid lines) have been corrected for 
source detection incompleteness. We also show these data without correction (dotted lines). In addition, we show (dashed line; black) the 
distribution of K-band optical light, and the best-fit simple powerlaw models to each dataset. In general, the GC distribution is flatter 
than the optical light, and the blue GCs tend to have a flatter distribution than the red GCs. 

This trend cannot simply be attributed to observational 
biases, such as the systematic failure to detect faint (low- 
mass), red (metal-rich) GCs, since that would actually 
result in a stronger correlation in the colour-magnitude 
diagram. Similarly, our assumptions of a constant age 
for all GCs or a blue horizontal branch for all GCs do 
not appear to bias our results. 

The apparent contradiction of no colour-magnitude 
correlation implying a mass-metallicity relation can be 
understood in terms of the weak dependence of the SSP 
mass-to-light ratio on the metallicity. To put this an- 
other way, if there was no correlation between mass and 
metallicity, due to the dependence of the mass-to-light 
ratio on metallicity, we would expect there to be a sig- 
nificant colour-magnitude correlation. To investigate this 
further, we simulated an artificial dataset comprising a 
similar number of GCs to that observed. For each ar- 
tificial GC, we randomly assigned a mass and metal- 
licity value from, respectively, uncorrelated log-normal 
and normal distributions which approximately matched 
the observed data. Using the appropriate V and I-band 



mass-to-light ratios from Maraston (2005), these data 
were converted to V and 1-band photometry, and esti- 
mated experimental noise was added. There was strong 
evidence of a weak correlation (the probability of no cor- 
relation was < 10~ 40 ) between colour and magnitude, 
with the fainter GCs being, on average redder. Con- 
versely, we have also used similar simulations to generate 
colour-magnitude diagrams from artificial GC datasets 
which obey the observed mass-metallicity dependence. 
Exactly as expected, we find that there is only weak ev- 
idence of a correlation in the resulting data (prob of no 
correlation ~0.1%, as compared to ~ 10~ 37 in colour- 
metallicity space), demonstrating that the metallicity- 
dependent mass-to-light ratio can, indeed, destroy a 
colour-magnitude relation. 

The galaxies in our sample span a range of magnitudes 
and so it is important to consider whether the known cor- 
relation between the galaxy mass and the peak of the GC 
colour distribution may affect our results. In part the re- 
lation reflects the increasing importance of re d clusters 
in higher-mass galaxies (e.g. Peng et al. 2006 ) and it is 
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difficult to imagine how that could introduce a spuri- 
ous GC mass-metallicity correlation if one is not a lready 
present. More problematically though, Peng ct al. report 
a correlation between the mean metallicity of the blue 
GCs and the galaxy mass and a similar trend for the red 
clusters. However, lower-mass galaxies are obse rved to 
have s ystem atically narrower and fainter GCLFs ( Jordan 
et al| |2007D , which would be consistent with a trend for 
brighter GCs to be redder. In any case, for the magni- 
tude range of the galaxies we consider here, the scatter 
about each correlation is at least as large as the trend. 

After this work was sub mitted f or pub lication, we be- 
came aware of a paper by Kundr ( 2008 ) which suggests 
that an apparent mass-metallicity relation may arise due 
to observational biases. He proposes that a weak mass- 
radius relationship, coupled with differing errors in the 
photometric aperture correction (since larger clusters 
may be marginally resolved) in different wavebands may 
lead to a spurious correlation. To investigate whether 
this could plausibly explain our observations, we first es- 
timated the error in the aperture correction for both V 
and I band by convolving King models with varying half 
light radii (R1/2) with PSF images (as described in § ||). 
Since the vast majority of the GCs in our study were 
detected on the WF CCDs and in galaxies more distant 
than 15 Mpc, we computed the correction for our WF 
photometry only and adopted a 15 Mpc distance while 
allowing R4/2 to vary from 1 to 8 pc (consistent with ob- 



relation is slightly less steep and has a slightly higher 
normalization. The observed trend for GCs may arise 
from a number of effects, such as self-enrichment (e.g. 



servations, e.g. Spitler ct al. 2006). Although the size of 
the cluster can affect the photometry in either band by as 
much as ~0.3 mag, the effect was very similar in each fil- 
ter so that V-I only changed by 0.1 mag. Such a small 
effect is not enough to wipe out the expected V-I versus 
V correlation in the case that mass and metallicity are 
uncorrelated. We have verified this by adding a colour- 
size relation into our Monte-Carlo simulations detailed 
above. We adopted a pathological monotonic linear re- 
lation between GC radius and V-band magnitude, which 
is much stronger than the weak trend observed. For 
our simulated colour-magnitude diagram under the as- 
sumption of no mass-metallicity relation, we found that, 
although the mass-radius effect slightly reduces the re- 
sulting colour-magnitude correlation, the probability the 
data are consistent with no correlation is still ~ 10~ 32 . 

While our results depend to some degree on the 
adopted SSP models , it is clear that t wo independ ent sets 
of m odels — those of |Marastorj| (2005 ) and those of Fioc & 
Rocc a- Volmerange ( |1997 ) both imply a significant corrc- 
lation between the mass and metallicity in these objects. 
The different normalizations and slopes for these models 
simply arise from differences in the predicted relations 
between the V and I-band mass-to-light ratios and the 
metallicity between the sets of models. However, any 
prediction of a metallicity-dependent mass-to-light ratio 
will clearly give rise to a mass-metallicity correlation of 
this kind. 

It is also interesting to compare our best-fitting mass- 
metallicity relation with that observed in early-type 
galaxies. We show in Fig || an extrapolation of the best- 
fitt ing straight- lin e relat ion for early-type galaxies found 



Mieske et al. 2006| ), in which the deeper potential wells 



of more massive GCs are more able to hold onto met- 
als ejected from the very first generations of stars within 
them than in less-massive systems. Alternatively, the 
mass of a forming GC may be related to t he mass (and 



meta llicity) of its parent gas cloud (e.g. [Harri s et al 
2006). Whatever processes actually underlie tli 



by Thomas et al. ( 2005|) . Intriguingly, we find a reason- 



able qualitative agreement, although the Thomas et al 



metallicity relation, the remarkable amount of intrinsic 
scatter (0.6 dex in metallicity) implies that stochastic 
processes actually dominate in determining the metallic- 
ity of an individual GC. Alternatively if there are, in fact, 
metal-rich and metal-poor sub-populations that cannot 
be distinguished in our data due to the statistical noise, 
and which exhibit different relations, the measured scat- 
ter will actually substantially overestimate the magni- 
tude of such effects. 

Recent (B,I)-band photometry, which allows better 
separation of metal-rich and metal-poor clusters, has re- 
vealed colour-magnitude correlations for the metal-poor 
(blue) sub-population of GCs in some early-type galax- 
ies, but no clear trend for the red population. Trans- 
forming a composite colour-magnitude diagram of three 
giant Virgo e l liptic als onto the mass-metallicity plane, 
Mieske et al. ( 2006 ) found little evidence of bimodality, 
but a trend of increasing average metallicity with mass, 
broadly consistent with our observed mass-metallicity 
dependence. Although bimodal colour distributions need 
not necessarily imp ly bimodal metallicity distributions 
(e.g. Richtlcr 2006| ), it is unclear whether the absence 
of bimodality in the mass-metallicity plane is simply an 
artefact of noise introduced during the transformation 
from colour-magnitude space. In particular, for our data, 
the error-bars on [Z/H] were typically so large that two 
putative distinct populations may have been blurred to- 
gether. This is particularly relevant when we consider 
whether the "blue tilt" and, in particular, the lack of 
a colour-magnitude relation for the metal-rich clusters 
found by previous authors, are consistent with our data. 
To assess this, we simulated a set of artificial GC pho- 
tometry data-points in the (g,z) photometric bands, with 
random z-band absolute magnitudes ranging from -7 to 
-11 and colours obeying the " blue tilt" relation found 
for blue GCs by Btrader et al. (2006). Transforming the 
data to the mas s-metallicity plane using the models of 
Maraston ( 2005 ), the data were approximately consis- 
tent with a simple linear relation of the form [Z/H]=- 
4.43+0.501og ln (M/M ( T ) ). Performing the same analysis 
for the red GCs found by Strader et al. , which we approx- 
imate as having g-z colours of 1.4, we found [Z/H]~-0.22. 
For a population of both blue and red GCs, we would 
expect the mean mass-metallicity relation to be simply 
a weighted average of these two relations. The existing 
best-fit slope and intercept are actually reproduced fairly 
well (0.25 and -2.3, respectively) for a (reasonable) blue 
GC fraction of ^50%. This excellent agreement with 
our results indicates that we may be observing the same 
trend as was observed in the (g,z) band, albeit the differ- 
ent trends for the red and blue populations are averaged 
together. 
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APPENDIX 

FITTING THE GCLF WITH A RESPONSE MATRIX 



Our method to fit the GCLF took account of both the source detection incomplete ness and th e effects of the 
Eddington bias by adapting the "response matrix" formalism from X-ray astronomy (e.g. Davis 2001). Wang] ( |2004 ) 



used a similar, but not identical, method to correct for the Eddington bias when fitting luminosity functions of X-ray 
binaries. For each galaxy we have a set of measured GC magnitudes (and intrinsic errors on the measured magnitudes), 
which we binned into an arbitrary set of bins. The resulting histogram can be written as Di (i = 1. . . n). Consider a 
GC with an intrinsic magnitude m, which we measure and assign to one of these bins. Due to statistical fluctuations 
there is a finite possibility that it will be detected in any one of the bins 1. . . n (or not at all). We can therefore define 
Pi(m) as the probability that this object is detected in bin i. If the differential GCLF can be written dn/dm, it follows 
that: 



f , i ^ din , 
< Di >= I dm Pi(m)— — (m) 

J-oo dm 



where the angle brackets denote the expectation. We can approximate this integral as a sum: 



< Di >~ ^Tp tj I dm -7— M = ^Tpijfj 



3=1 3 1 j=l 

where pij = pi(n%j), and the spacing between each rrij (j = 1. . . N) is small enough to ensure that Pi(m) does not 
vary strongly between them. For any given parameterized GCLF, w e therefore have the expected number of GCs in 
each bin i=l. . . n, allowing us to construct a likelihood function after |Cash (1979). Maximizing the likelihood function 



enables the parameters which define the GCLF shape to be constrained by the data 

The probability p^ is, in principle, dependent on the measurement errors, source confusion and the "background" 
level (which varies strongly across the image) so its computation from first principles is extremely demanding. Instead, 
we choose to estimate it by adding objects of known intrinsic magnitude mt to the image and measuring them from the 
data. Since this is equivalent to setting fj = djk, where Sjk is the Kronecker delta function, it fo llow s that < Di >= p^. 



For each magnitude rrif., we added 100 stars in this manner to the image (as described in § 4.1), and assumed that 
the arithmetic mean of our measured Di recovered for those stars was a good approximation for the expectation in 
order to estimate p^.. The accuracy of this assumption is dependent on the number of simulations performed and so 
we expect some intrinsic inaccuracy in our measured probabilities which, in turn, may affect the likelihood function 
and derived parameters. It is therefore important to assess whether 100 simulations are sufficient for our purposes. 

To assess the accuracy of the matrices pij, we performed a set of Monte-Carlo simulations to determined how 
accurately a known GCLF can be recovered. We adopted an arbitrary analytical expression for the incompleteness as 
a function of measured magnitude, which was chosen approximately to match the behaviour found for a representative 
galaxy in our sample, NGC4365. We stress that we do not have to match exactly the properties of the galaxy, only 
approximately for our present purposes. Adopting measurement errors consistent with the observed data, we were 
able to simulate an "observed" dataset, given a known input luminosity function. In this way, we were also able to 



generate corresponding response matrices, follo wing the procedure described in § 4.1. We generated matrices using 



100 simulated sources per magnitude bin, as in § 4.1 but also generated a separate set of matrices using 10000 sources 
per magnitude bin. We subsequently refer to these as the standard and large matrices, respectively. 

We carried out 100 Monte Carlo simulations, on each of which we generated both a dataset and corresponding 
response matrices, binned identically to the real data. For the input GCLF, we assumed my= 24.0, ay = 1.3, 
consistent with NGC4365, and a total number of 500 GCs. For each dataset and response matrix pair, we fitted the 
simulated data, allowing both niy and the normalization to fit freely. The mean and standard error of the best-fitting 
parameters, averaged over the simulations, can be used to assess the bias and statistical error (which will also include 
a contribution from the error on the response matrix due to the finite number of simulated sources used to create it) 
implicit in our analysis. Applying the standard matrices, we obtained best-fitting niy= 24.00 ± 0.12 (1-cr error) and a 
normalization of 504 ±41 (1-cr) GCs. The best-fitting values are sufficiently close to the true value to indicate that the 
results are not biased, and the recovered errors are broadly consistent with those found for the real data. Adopting 
the large matrices, we recovered best-fitting values of my= 24.02 ± 0.11 (1-cr errors) and a normalization of 506 ± 36 
(1-cr) GCs. Since increasing the number of simulations used to generate the matrices by a factor 100 has virtually no 
impact on the recovered parameters or error bars, it follows that the standard matrices are sufficient to fit the GCLF. 

As another way to assess the accuracy of the matrices, we also performed a set of simulations for each dataset, 
in which the response matrix was replaced with an artificial matrix, p'^, which was derived from the real matrix, 
with approximately the expected level of statistical noise added (assuming the "observed" p^ is exactly correct), and 
the data refitted. We performed 25 such simulations for each object and adopted the 1-cr scatter of the best-fitting 
parameters from our Monte-Carlo simulations as an estimate of the uncertainty arising from our matrix calculation. 
We computed p'^ in the following manner. For each magnitude rrij, we drew 100 sources which were distributed 
amongst the histogram data-bins i=l...n randomly but following the probability distribution p^. Denoting the 
resulting data histogram as (D'^j, clearly < (D'^j >= 100pij (since there are 100 sources). The actual values of (-D-)j 
will be distributed about this expectation due to statistical scatter, allowing us to define p'^ = (D'^j/IOQ. In practice, 

we found the uncertainty arising from the matrix computation for Sn, Sl and my to be < 10%, <4% and < 1%, 
respectively, which are far smaller than any other source of error in our calculations for these quantities. 



14 



Humphrey 



TABLE 5 

Catalogue of GC candidates 



Source 


RA 


dec 


V 


I 


Mass [Z/H] 








mag 


mag 


10 6 M Q 



NGC 1332 

1 03h26ml7.34s -21°20'25.5" 21.93 ±0.04 20.79 ±0.03 

2 03h26ml6.94s -21°20'16.7" 23.9 ± 0.2 22.9 ± 0.2 

3 03h26ml7.38s -21°20'18.4" 23.4 ±0.1 21.99 ±0.09 1.11 ±0.07 

4 03h26ml7.58s -21°20'18.9" 22.64 ± 0.08 21.36 ±0.06 2.08±o2 9 
5 03h26ml7.82s -21°20 / 19.1" 23.4 ± 0.1 22.22 ± 0.09 0.69 ±0.15 -0.22 



2 6+ ' 5 
n 9 o+0.27 



-0 23+ 15 
i O+0.6 

fi7 +o.oio 

u -°' -0.29 
44+0' 2 3 

+0.77 

-n.sn 



Note. — The complete version of this table is in the electronic edition of the Journal. The 
printed edition contains only a sample. 

SOURCE LISTS AND PHOTOMETRY 

We show in Table [s] a list ol the GC candidates detected in each galaxy. For each source, we provide a unique 
identification number within that galaxy, as well as its right ascension and declination. We provide the V and I-band 
photomet r y for e ach candidate GC as well as the mass and metallicity, having converted them using the SSP models of 
|Maraston| (|2005|) , assuming an age of 13 Gyr and a blue horizontal branch. All errors are quoted at the 1-er confidence 
level. Note that [Z/H] is constrained to lie in the range covered by the models, i.e. -2.25 to 0.67. 
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